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Abstract. The momentum distributions associated with both the wave function 
of a particle behind a grating and the corresponding Bohmian trajectories are 
investigated and compared. Near the grating, it is observed that the former does 
not depend on the distance from the grating, while the latter changes with this 
distance. However, as one moves further apart from the grating, in the far field, 
both distributions become identical. 

1 Introduction 

Quantum interference experiments where beams of one particle at a time are used have intensi- 
fied the theoretical investigation of the topology of particle trajectories behind the interference 
grating 1 2345678 . The aim of all the approaches developed in this direction is to get con- 
sistency between the quantum mechanical particle distribution and the distribution associated 
with the trajectories of the particles. In this paper we compare the features of Bohmian trajecto- 
ries with those displayed by the trajectories determined using the momentum distribution (MD 
trajectories) associated with the wave function of the particle. As is well known, Bohmian tra- 
jectories follow the flux lines of the quantum flow, and therefore reproduce exactly the quantum 
mechanical distribution in both the far and near field 1 2 3 4 6,8]. In particular, it is remarkable 
the consistency of a set of Bohmian trajectories in the near field behind a multiple slit grating 
within the context of the Talbot effect [8] . This effect has been observed experimentally with 
relatively heavy particles, such as Na atoms [9] or Bose-Einstein condensates |10j . 

Moreover, here we also study the momentum distribution associated with Bohmian tra- 
jectories, and compare it with the momentum distribution associated with the wave function 
of the particle. We find that in the far field these two distributions are identical. However, 
in the near field, the distribution of transverse momenta associated with the Bohmian trajec- 
tories changes with the distance from the grating, and therefore differs from the momentum 
distribution associated with the wave function. 

The essential feature of the Bohmian deterministic trajectories is that a particle passing 
through different slits will never reach the same point on the detection screen 0]. On the 
contrary, this limitation does not hold for MD trajectories arising from different slits, which 
may reach the same point on the screen |5|7) . The momenta of particles moving along MD 
trajectories are distributed according to the momentum distribution determined by the wave 
function of the particle. The probability of arrival of a particle to a given point on the screen 
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depends on: (1) the slit through which it passed and (2) the presence of other slits. The second 
type of dependence is a consequence of the dependence of this probability on the transverse 
momentum distribution, which depends on the properties of the grating as a whole [5]. As 
shown, MD trajectories reproduce fairly well the quantum mechanical space distribution in the 
far field. 



2 Wave function of a particle behind a grating 

The motion of a particle behind a grating is determined by its wave function. Here, we assume 
that in front of the grating (y < 0) we have a plane wave with initial momentum p = Kk = hki y 
(and de Broglie wavelength A = 2ir/k), moving along the longitudinal direction y. Behind the 
grating (y > 0), the solution can be expressed as a product of the longitudinal part, a plane 
wave, and the transverse part, i.e., 

&(x,y,t)=Be ik y e - i " t TP(x,t), (1) 

where B is a normalization constant. Arsenovic et al. have shown [5] that the transverse part 
(for simplicity taken to be one-dimensional) is given by 

i r+oo 

$[x,t) = -=l dk x c(k x )e^ x e' ak ^ 2m , (2) 



where 



c(k x ) = -= I dx^(x,0)e~ l ^ x (3) 
V 2n 



is the Fourier transform of the initial transverse wave function tp(x, 0), which is the wave function 
just behind the grating, is related to the transmission or window function associated with the 
grating, and refers to the time of passage of a particle through the grating. Substituting ([3]) 
into (|2|), ip(x,i) is expressed in terms of the initial wave function as 



^,i) = 4=e^/ 4 ^{x'^y^-^l^dx 1 . (4a) 
v2Trht J -co 

By assuming that the motion along the y-axis is classical, i.e., it satisfies the relation y = vt = 
(hk/m)t, one finds the dependence of t()(x,t) on the longitudinal coordinate y: 

ip(x, t = ym/Tik) = -j= e-" /4 / ip(x', 0)e lk ^- x '^^ y dx' . (4b) 



/27n/ 

Expressions (|4"aj) and (|4b[) are particularly useful when ^(^,0) consists of discrete pieces of 
elements where it is zero. In the case of a one-dimensional grating with n slits of equal width 
S (see Fig. [IJ, from the boundary and normalization conditions follows 

^(au,0) = 



at the openings, and 

V>(x,o) = o 

elsewhere (i.e., outside the openings). Note that this is an idealization of the effect that a perfect 
periodic grating with full transmission in the openings would have on a plane wavefront that 
would reach it. Since all the points on this wavefront have the same phase, the transmitted 
pieces of wave function described above will also have the same phase. We can neglect this 
phase, since introducing a constant phase factor will not alter the physics of the problem. The 
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Fig. 1. Scheme of the physical sys- 
tem studied. Note that x' refers to 
the coordinate along the grating, 
while x is the actual coordinate out 
of the grating. 



probability amplitude of the particle transverse momentum c(p x ) — c(k x )/Vh — c(p x /h)/ \ffi 
is then given by 



c(k x ) = 



2 sin (Y) sin 
mS k x ginf^) 



(5) 



where d is the distance between the center of two consecutive slits. 

In Optics it is common to consider sharp-edge slits to study interference and diffraction, 
as in the model described above. Changing the shape of the initial wave function would not 
induce significant changes in the momentum distribution. For example, one may use Gaussians 
instead of square wave functions, i.e., at the openings we have 



0) 



~(x— x c ) 2 /a 2 



(6) 



where x c is the center of the corresponding opening and a is of the order of the slit width. Then 
a different multiplicative Gaussian factor appears in the first contribution to the momentum 
distribution, thus c(k x ) becoming 



c{k x ) oc e 



,2 2 , . sin 

— k x a /4 



/ k x dn \ 



sin 



(7) 



For x 3> x', the term x' can be neglected in the exponential function inside the integral in 
dH). By doing that, one obtains [5] 



TTl —iir/4 i x '" 

*> ' e 2nt 



2nht 



dx'^(x',0)e- 



(8) 



The integral in © is nothing else but the function c defined in ([3]), taking into account that 
the variable xm/ht plays the role of k x . Therefore, one can write 



rn 
tit 



—in/4 ix m/2ht 

e 1 e 1 C 



(9) 



The wave function |T]), where the transverse part is given by @ or ([?]), explains in a unified way 
many effects and properties of particle diffraction and interference [7] , in particular the Talbot 
effect. Note that, for a periodic (infinite) grating, the integral in |2]) reduces to the sum over 
discrete momentum values and one can easily see that at integer multiples of 2Lt (Lt = d 2 /X 
being the Talbot distance 8 9 ID]) the transverse wave function is simply equal to the wave 
function on the grating. The grating pattern is also repeated at odd multiples of Lt, but shifted 
by half a grating period, d/2, along x-axis. 

In the case of a finite grating, an approximate form of the Talbot effect appears [7111] , 
Namely, in the near field one finds self-images distorted at the ends. 
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3 Momentum distribution trajectories 

Interference experiments with beams such that only one particle is launched at a time have 
demonstrated that particles accumulate with time on the detection screen, building up an 
interference pattern after some time. In order to describe the emergence of the interference 
pattern through an accumulation of single particle events one has to assume that particles 
move along certain trajectories. 

The theoretical investigation of the topology of the trajectories is based on various assump- 
tions. Arsenovic et al. ,5] assumed that trajectories were straight lines starting at different 
positions on the slits. The longitudinal momentum of a particle moving along such a trajectory 
is equal to the longitudinal momentum of the incident beam that reaches the grating, and the 
distribution of the transverse momentum of the particles is given by |c(p K )| 2 . This is why we 
call these trajectories the momentum distribution (MD) trajectories. 

Within the MD approach, a particle can start from any slit; if it has the right transverse 
momentum, it will reach the chosen detection spot. The expression for the screen arrival prob- 
ability based on the idea of MD trajectories [5] reads as 

/+oo r + co 

dk x I dx'\c(k x )\ 2 \tp(x',0)\ 2 S(x- x' - hk x t/m). (10) 
-oo J — OO 

The total probability can be expressed as a sum of n terms 

n 

P(x,t)=^2Pi(x,t), (11) 

8=1 

where 

Pi{x,t) = — / \c{k x )\ 2 dk x (12) 

nd Jm{x-x%)/ht 

is the probability that a particle which passed through the i-th slit of the grating arrives to a 
point (x,y = Tikt/m) at time t. Here x) and x\ are the coordinates of the left and right edges 
of the i-th slit. Note that, as said in the previous Section, the momentum distribution c(k x ) 
depends on the initial wave function chosen. This will only influence the probabilities along the 
x direction, but not the shape of the MD trajectories. 

Numerical calculations for various gratings show [5111112) that the distributions P(x, t) 
and \4>(x,t)\ 2 are almost identical in the far field. In the near field, both distributions look 
qualitatively similar, however they differ numerically as well as in certain important details. 
This is understandable because near the slits the topology of particle trajectories is more 
complicated than in the far field, where the straight-line approximation works fairly well. 



4 Bohmian trajectories 



The Bohmian trajectories associated with the state of a single particle, \P(r,t), are determined 
PQ from the differential equation (guidance condition) 



dr 1 
v(r,t) = - = -V5(r,t), 
dt m 



(13) 



where S(r,t) is the phase of the wave function written in the polar form, i.e., 

W{r,t) = Re iS/n . (14) 

It can be easily shown that 

WS = Mm — . (15) 



5 



4.5 




-0.4 -0.2 0.2 0.4 

x[nm] 



0.002 
0.001 








: ! 








M 


I; \ ;}, : 


■;, f 








i 1 1 

11 I 

!| r| : 


/J. j 






4 


Ifi . M 
£ 1 


: JH -7: g; ;j< 






\ I, 


1 i ^ 
| ' !! ;: j 


;:; .; ■ .■ . ; ■ 
: ■ . .,: .. 























-0.4 -0.2 0.2 0.4 

X[um] 



Fig. 2. a) Particle distribution function be- 
hind a grating with n=5 slits for y—1.25Lr, 
where Lt is the Talbot distance Lr=d 2 /A. b) 
Trajectories behind a grating with n=5 slits 
for y < 1.251/t. Trajectories passing through 
different openings are represented with differ- 
ent line styles. The values of the parameters 
used in the calculations are: d = O.lxlO -6 m, 
8 = 0.05xf(T 6 m, m = 1.19xl0 -24 kg, 
v = 220 m/s and A = 2.53xlCT 12 m. 



Hence, according to JT]), in our case we have 

dr Ti 
dt m 



i y k + i x Im 



d x ip(x,t) 



where d x = d/dx. The equation of motion along the y-axis, 

hk 



is very simple, and has the solution 



dy 
dt 



hk 

y = yo H — t. 



On the other hand, the equation of motion along the :r-axis 



dx h 
—— = — Im 

dt m 



d x ip(x,t) 



ip(x,t) 



(16) 

(17) 
(18) 

(19) 



can not be solved analytically for all values of x. Unlike MD trajectories, for Bohmian trajecto- 
ries both their distribution and shape are influenced by the choice of the initial wave function 
[3], although fundamental features still remain. 

Equation (|19|) can only be solved analytically in the far field, using the approximate expres- 
sion ((9]) for ip(x,t). Using this approximation, Eq. (fT9]l reduces to 



dx 
~dt 



x 

7' 



whose solution reads 



to 



(20) 
(21) 
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Fig. 3. a) Particle distribution function be- 
hind a grating with n=5 slits for y—12.5Lr , 
where Lt is the Talbot distance Lr=d 2 /A. b) 
Trajectories behind a grating with n=5 slits 
for y < 12.5Lt. The values of the parameters 
used in the calculations are: d = O.lxlO -6 m, 
8 = 0.05xl(T 6 m, m = 1.19xl0~ 24 kg, 
* m l v = 220 m/s and A = 2.53x 1(T 12 m. 




We have found numerical solutions of Eq. (p~9|) and plotted the Bohmian trajectories for particles 
behind gratings with n = 5 and n = 30 slits. As can be seen in Figs. 1 and 2, the density of 
Bohmian trajectories at a certain distance from the grating is in good agreement with the 
quantum mechanical probability density, \ip(x,t — ym/hk)\ 2 , in both the far field (see Fig. 2) 
and the near field (see Fig. 1). Bohmian trajectories, grouped in bunches, follow the directions 
that end at the regions of the different intensity peaks. Such an agreement was found previously 
by Sanz et al. [11314] . who plotted together the intensity pattern obtained by means of the 
standard quantum mechanics and the histogram obtained by counting Bohmian trajectories 
(see Fig. 5 in 0], for instance). 

The consistency of the set of Bohmian trajectories in the near field within the context of 
the Talbot effect found by Sanz and Miret-Artes [H] is remarkable. This consistency is also seen 
in Fig. 3, which shows the Bohmian trajectories behind a grating with n = 30 slits (only one 
half of the full pattern is shown; the other half is its mirror image). 



5 Comparison between Bohmian trajectories and MD trajectories 



An essential feature of the Bohmian deterministic trajectories is that a particle passing through 
different slits will not reach the same point on the detection screen. In the far field, Bohmian 
trajectories asymptotically approach straight lines that connect the center of the grating with 
the detection spot (see Fig. 4). On the contrary, MD trajectories from different slits may reach 
the same point on the screen, because this trajectories are associate with different values of 
transverse momentum. Namely, a particle which leaves certain point at the slit may have various 
values of momentum, in accordance with the transverse momentum distribution. 
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Fig. 4. Bohmian trajectories behind one half of a Ronchi grating with n=30 slits. The values of the 
parameters used in the calculations are: d = 0.2xl0~ 6 m, S — O.lxlCP 6 m, k — (7r/8)xl0 12 m -1 , 
m = 3.8189xl0" 26 kg and v = 1084 m/s. 



Fig. 5. Sketch of MD (five straight lines) and Bohmian 
(curve) trajectories where there essential features are 
compared. The central among these straight lines is the 
asymptote of a Bohmian trajectory. 
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Fig. 6. Histogram of the Bohmian momentum distribution at four distances from a grating with 
n=5 slits: (a) j/=Lt/40, (b) y—Lr/i, (c) y—1.25Lr and (d) y—12.5Lr ■ In panel (d), the quantum 
momentum distribution has also been plotted with dashed line. 



6 Momentum distribution 



It is evident that the probability amplitude of the transverse momenta is an essential feature 
of the wave function and wave field behind a grating. Thus, the following question arises: is the 
distribution of transverse momenta associated with Bohmian trajectories identical or different 
from the distribution |c(p K )| 2 = \c(k x )\ 2 /ft,, where c(k x ) is determined by J^J? 

Long time ago, based on a general expression for the distribution of Bohmian momenta, Tak- 
abayasi [13] concluded that the aforementioned two distributions were different functions. We 
have found that this conclusion is not always true. Using the Eq. (|19p for Bohmian trajectories, 
and the approximation ([9]) for the wave function, we are going to show that the distributions 
determined from the wave function and from Bohmian trajectories are identical in the far field. 
In the near field, the distribution of transverse momenta associated with Bohmian trajectories 
changes with the distance from the grating and is different from the distribution |c(pa;)| 2 . 

An analytical expression for the transverse momentum distribution in the far field can be 
easily obtained. One starts from a general relation between the probability distribution in the 
x-space and in the p^-space, 

P(x,t)dx = n(p x ,t)dp x . (22) 

Taking into account that P(x,t) = \ip(x,t)\ 2 as well as relation (|19p . we find the following 
expression for the probability distribution of Bohmian momenta along the x-axis: 

nM = -teT = -fcs-> (23) 
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Fig. 7. Histogram of the Bohmian momentum distribution at four distances from a grating with n—30 
slits: (a) y=L T /i, (b) y=L T /2, (c) y=0.75L T and (d) y=L T . 



where d xx = d 2 /dx 2 . Using the approximate expression © for the transverse wave function, 
which is valid in the far field, we find: 

d 2 S _ m 
dx 2 t 

Therefore, the probability distribution of Bohmian momenta in the far field will be 



(24) 
(25) 



n{p x ,t) = ^\c(k x )\ 2 = \c(p x )\ 2 . 



(26) 



We do not have an analytic expression for II(p x ,t) in the near field, but the histograms of the 
probability density of Bohmian momenta plotted in Figs. 5 and 6 show that it changes with 
the distance from the grating. This means that, in the near field, the probability density of 
Bohmian momenta differs from the momentum distribution |c(p x )| . 



7 Conclusions 



By appealing to descriptions in terms of quantum particle trajectories in interference experi- 
ments, one can understand how interference patterns emerge from the accumulation of single 
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particle events. This justifies theoretical investigations of properties and statistics of particle 
trajectories. 

In Bohmian mechanics trajectories reproduce exactly the quantum mechanical space distri- 
bution in both the far and near fields. The consistency of the set of Bohmian trajectories in the 
near field within the context of the Talbot effect, demonstrated by Sanz and Miret-Artes [5], is 
remarkable. 

The probability amplitude of transverse momenta is also an essential feature of the wave 
function and the wave field behind the grating. Taking this into account, the following ques- 
tion arises: should particle trajectories comply with the transverse momentum distribution? 
From our numerical calculations and analytical treatments, it follows that the distributions of 
transverse momentum determined from the wave function and from Bohmian trajectories are 
identical in the far field. On the other hand, in the near field, the distribution of transverse 
momenta associated with the Bohmian trajectories changes with the distance from the grating 
and is different from the distribution |c(p x )| 2 . 

Considering that the answer to the aforementioned question should be positive, Arsenovic et 
al. proposed [5] to approximate trajectories by straight lines and to assume that the distribution 
of particle momenta is determined by the wave function. These trajectories, which we call MD 
trajectories, reproduce fairly well the quantum mechanical space distribution in the far field 
[5111] , However, in the near field the agreement is not so satisfactory. It seems that a better 
agreement of the space distribution derived from MD trajectories with the quantum mechanical 
space distribution in the near field could be obtained by combining peaces of various Bohmian 
trajectories |H2l3l4l8j . and by studying also lines of a quantum mechanical current. 
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